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Abstract 



> 

\^q \ The normalized factorial moments F q are continued to noninteger val- 

ues of the order q, satisfying the condition that the statistical fluctuations 
remain filtered out. That is, for Poisson distribution F q = 1 for all q. The 
continuation procedure is designed with phenomenology and data analysis 
in mind. Examples are given to show how F q can be obtained for positive 
and negative values of q. With q being continuous, multifractal analysis 
is made possible for multiplicity distributions that arise from self-similar 
dynamics. A step-by-step procedure of the method is summarized in the 
O . conclusion. 

> 

•i-H 

1 Introduction 

Continuation of the normalized factorial moments F q to arbitrary, noninteger values 
of q is not just a mathematical problem. It has high phenomeno logical significance, 
and provides a powerful method to analyze experimental data that can reveal aspects 
about multiplicity fluctuations hitherto unexplored. 

A historical review of the problem is in order. Bialas and Peschanski [0 first 
introduced F q as a means of studying scaling behavior of multiplicity fluctuations as 
a function of the resolution scale 5, a subject usually referred to as intermittency. 
Although a large part of the effect seen in the data has recently been found to be 
due to Bose-Einstein correlation among like-charge particles 0, [3], intermittency at 
a weaker level is still present for the unlike-charge particles and its origin remains to 
be clarified ||, especially for q > 2. There is far more dynamical information about 
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high-energy collisions than can be uncovered by studying two-particle correlation only. 
With that point of view forming the basis of our discussion here, we now outline the 
problems associated with the use of F q . 
Let us first recall the definition of F q : 

(n{n-l)---{n-q+l)) 

F « = w {) 

where (■ ■ •) denotes an average weighted by the multiplicity distribution P n . The most 
outstanding property of F q discovered in Ref. Jj]] is that it filters out the statistical 
fluctuations, so any nontrivial behavior of F q is a direct indication of some features 
about the dynamics of particle production. A quick review of that will be given in the 
beginning of the next section. Another significant aspect about F q is that an event 
can contribute to ([I]) only if n > q. Thus for small 5 where (n) is small in a bin, only 
rare events with high spikes (n > q) contribute. That is why it is sometimes said that 
intermittency measures spiky events. There are, however, disadvantages in the use 
of F g . A corollary to the ability to select spiky events is its inability to extract any 
dynamical information about dips. It is by now generally recognized that rapidity 
gaps, like voids in galactic structure, are important to study. Those are, of course, 
large dips. In nuclear collisions where multiplicity per bin is large, unusual dips, which 
can be small but deep, are as significant as unusual spikes. For such fluctuations it is 
necessary to study F q for q < 1, especially negative q. Furthermore, for multifractal 
analysis of multiparticle production the continuation of F q to noninteger values of q is 
necessary in order to allow differentiation with respective to q. These studies cannot 
be done, if F q is defined as in ([!]). 

A method to investigate moments of arbitrary order was suggested several years 
ago in terms of the G moments ||. It was later modified to achieve better power- law 
behavior ||. However, in overcoming the defects of F q , the G moments fail to retain 
the principal attribute of F q , i.e. the screening of statistical fluctuations. Subtraction 
of the statistical component has to be done by hand . Since that can be achieved 
only by simulation, the method is not elegant. But it has provided the first glimpses 
into the multifractal structure of particle production. 

In this paper we describe a method that retains both attributes: it eliminates 
statistical fluctuations and is defined for continuous q. Although the mathematical 
technicalities involved may at first sight appear to be of theoretical interest only, the 
method is developed with phenomenology in mind. The purpose of the program is to 
extract quantitative information about dynamical fluctuations from the experimental 
data and to present it in a form suitable for comparison with theoretical predictions. 
Thus the problem of data analysis has not been overlooked in favor of mathemat- 
ical expediency in the hope that the procedure can be readily accessible to direct 
experimental application. 



2 The Problem 

We first review the virtue of factorial moments. To say that F q filters out statistical 
fluctuation, one first assumes that the latter enters the multiplicity distribution as a 
convolution with the dynamical component 

P n = S®D (2) 

where S represents the statistical component, which we take to be the Poisson distri- 
bution V^\ and D is the dynamical distribution. More specifically, (0) implies 

poo f n 

P n = / dt-e-'Dit) , (3) 

jo n! 

Let the numerator of (jl]) be denoted by f q , i.e. 

'•=£(^>f p - • (4) 

which is well defined for q being a positive integer. Substituting @ in (]D and 
performing the summation yield 

poo 

f q = / dtt«D(t) . (5) 



This is the gth moment of the dynamical D(t) and is free of statistical contamination. 
Since f\ = (n), we have 

F q = u/n . (6) 

If P n is a Poisson distribution, then D(t) = S(t — (n)) and f q = (n) q so 

F q = 1 , for integer q > 1 . (7) 

Because of this trivial result, one can state that any nontrivial F q reveals the nontrivial 
properties of D(t). 

To generalize F q to noninteger q, it must first be recognized that there is no unique 
continuation to complex q. Since P n must vanish as n —> oo (in fact, it must vanish 
for n > N for some finite iV at any finite collision energy), there are only a finite 
number of the F q moments defined at integer q values. Without an accumulation of 
F q at infinite q, unique continuation to noninteger q is not possible. Put differently, 
we can add to F q any arbitrary function that vanishes at the finite range of integer q 
where F q is specified and generate another function F q at noninteger q. 

A simple way of continuing (^) to arbitrary q is to replace the factorials by gamma 
functions, i.e. 

^o r(n-g+l) 



Procedures similar to (|5p, such as continuing d q G(z)/dz q to fractional q M, have been 
considered previously ||-[|Il[]. Since V(z) has poles and oscillates rapidly among those 
poles when z < 0, F q as defined in @ oscillates between large positive integers of q 
and is highly suppressed at large negative q. The question is whether one wants that 
kind of behavior at noninteger values of q. If not, what are the guidelines by which 
one makes alternative choices of the continuation schemes? 

In our view the only guideline is the primary rationale for considering factorial 
moments in the first place. And that is the elimination of statistical fluctuation at 
all q, not just at integer values of q. If one substitutes Poisson distribution into (||), 
one will find that F q ^ (n) q . Fig. 1 (a) and (b) show the results for F q when 

P„ = Pi°> = ^Pe-<»> , (9) 

for (n) = 6. Although (0) remains true for positive integers of q, F q is by no means 
equal to one for all q. The oscillations have larger amplitudes at high q, as revealed 
in Fig. 1(b). The situation is worse at smaller (n). In Fig. 2 we show the result for 
(n) = 1 in (|9]), for which (|) and (§) are again used. Notice that there is no connected 
region of q in which F q = 1, not even between q = 1 and 2. Between q = 9 and 10 
the peak of \F q \ is greater than 10 4 . From these results it is therefore not possible 
to claim that (f|) contains no statistical contribution at noninteger q. The reason for 
studying F q at noninteger q is consequently lost. The G q moments []5], would be 
better. 



3 The Solution 

To achieve our aim of retaining only the dynamical fluctuation in our continuation 
procedure, we demand that (|5|) defines f(q) for all q. Hereafter, we use the notation 
that when q appears as an argument of a function, instead of as a subscript, it is to 
be regarded as a continuous (complex) variable. A consequence of that condition is 
that the normalized factorial moment function satisfies 

F(q) = 1 (10) 

for all q in the case of Poisson distribution. Equation ( |10"D should be used as a test 
of the continuation procedure. 

The burden of this procedure is to determine D(t). If a theory specifies the 
dynamical distribution D(t) completely, then (|5|) prescribes a unique continuation 
of f q to the complex function f(q) at any q. But how is that to be checked by 
experiments where only P n is measured? Thus the procedure must supplement ([5|) 
with a way of determining D(t) from P n . This deconvolution process also cannot be 
made unique. The discrepancies show up as deviations of F(q) from 1. The region 
where ([TUP fails significantly can fortunately be controlled and pushed to large \q\. 



To deconvolute (|3]) one could consider making the inverse Laplace transform of the 
generating function G(z). However, there are difficulties connected with the fact that 
G(z) determined from the experimental P n is a polynomial having no singularities in 
the finite z plane. 

Our proposal is to expand P n in terms of negative binomial distributions (NBD) 
P^ B (j). One of the attributes of NBD is that it can also be expressed as a Poisson 
transform ||12| , as in (||]). They do not form a complete set of orthogonal functions, 



so in general they cannot be the basis functions for the expansion of an arbitrary 
function. However, we do not have an arbitrary function. The experimental P n 
(ignoring errors for the moment) is a set of N + 1 numbers for n = 0, 1, • • • , N. Thus 
the expansion 

^ = E^f(j) (ii) 

3=0 

is well defined with N + 1 coefficients a,, provided we specify P^ B (j) appropriately. 
One could consider other distributions instead of NBD, but for factorial moments 
that we shall eventually calculate NBD is most convenient. 
Now P^ B (j) is defined by p 



■)NB (i \ i. (fl -\- Kj) I Kj \ I Xj 



e?" &,*,) = r r iJ-'A t^- t^- , (12) 



T(n + l)T(kj) \kj + Xj J \kj 



x 



j , 



where Xj and kj specify the mean and inverse width of P^ B {j). How they depend on 
j will be discussed in the next section. We remark that Xj is equal to 

oo 

n(j) = T,^Pn B U) (13) 

n=0 

only if the sum extends to oo. Since our method is designed with phenomeno logical 
analysis in mind, where P n is given only for n — 0, • - • , N, all sums over n will be from 
to N, whether the summand involves P n or P^ B ■ Consequently, Xj is not exactly 
n(j). Although the discrepancy is small for the Xj and kj to be chosen, accuracy will 
be important, as we shall see. Extending the sum in (|TTD to a larger upper limit N' 
with P n = for N + 1 < n < N f would cause ctj to be very large and highly sensitive 
to the accuracy of the calculation; it is a procedure that should be avoided. Hereafter 
N will always be the maximum value of n for which P n is measured to be nonzero, 
and Xj and kj should only be regarded as real parameters of P^ B (j) that will be 
varied in the expansion in (11). It should be recognized that because P^ B (j) are all 
positive (unlike a harmonic function) and small on the wings, cij will have alternating 
signs and can have large absolute values if P n is small for a range of large n values. 
Thus accuracy in the ensuing calculations will be essential. With Pn B (j) specified, 
(|TTD is a set of iV + 1 simultaneous algebraic equations that can be solved for Oj in 
terms of the experimental P n . 



Define D NB (t) by the negative-binomial versions of (0), i.e., 

roc t n 



P n NB U)= rdtLe-Wfrj). (14) 

JO TV. 

Then it is known that \T2\ 



kj j-ki — l 



^^-(fj'm^'"- (15) 



The substitution of (y) and (|Tj) to the two sides of (|TTj) results in 

N 

D(t) = Y,a 3 D NB (t,j). (16) 

3=0 

Thus we have inverted @ and extracted the dynamical distribution D(t) from the 
experimental data on P n . In principle, this D(t) can be compared directly with the 
theoretical distribution. However, the more familiar arena for comparison involves 
the factorial moments, which are more closely related to the data. 

An interesting side remark that can be made here is that the feasibility of deter- 
mining D(t) from the data makes possible an experimental look at such theoretical 
quantity as the Landau free energy if the data on hadronic multiplicity distribution 
correspond to quark-hadron phase transition |13|| , or if the data are on photon dis- 
tribution at the threshold of lasing in quantum optics |IJJ. That is because in such 



problems D{t) in (|3]) is e~ F ^\ where F[t] is the free energy of the system. 

Returning to the problem on the factorial moments, we substitute ([R]) into our 
basic equation, (|^), for continuation to complex q, and get 

N 

M = T / ^f NB ^J) (i7) 

3=0 

where 

'"<'■•« - (I)" ' rJ wr ■ < 18 > 

This expression is obtained by integration over t in (^) and is valid only for 

Req>-kj. (19) 

Thus the domain of q that can be continued into before encountering the first singu- 
larity is governed by the smallest value of kj. 

It should be noted that whereas f q as defined in (f|) is obtained directly from the 
input P n , f(q) is determined from ( |17|) after the inversion is done and the continuation 
procedure followed. Even at positive integer values of q, f(q) is not exactly equal f q 



because of finite accuracy in the computation. Following (||), we define the continued 
normalized factorial moments by 



F(q) = f(q)/f(iy, 
instead of f(q)/(n) q . From (|T7|) and ( |18D we then have 



N 



m = E 



3=0 






where 



F™(q,j) 



/(I 

_ r(g + kj 
T(kj)kj 



F NB (q,j] 



(20) 



(21) 



(22) 



Equation (21) is our result, which is explicit once the values of aj are determined. 



4 The Range of x 2 and kj 

To complete the description of our continuation procedure, it is necessary to specify 
Xj and kj, which represent the average and inverse width of P^ B (j). They should be 
chosen to be not too far from the values of x and k of the input P n , where 



N 



x = (n) = ^J n P n 

n=0 



(23) 



k = (F 2 - 1) 
We therefore define first 

A, 



{n{n — 1)) jx 2 . 



Al-l + l 
1 2 N 



(24) 



(25) 



which ranges from —A/2 to +A/2 in equal steps, as j varies from to N. For Xj 
and kj to vary from the lower to higher sides of x and k, we set 






x(l + Aj) , 



(26) 



kj = k{l + Aj). 



(27) 



Thus Pn B {j) with lower Xj has wider width (lower kj); while that with higher Xj has 
narrower width. If A is not too large, all components Pn B {j) will have comparable 
magnitudes at large n, which should be small where P n is small. Otherwise, a,j 
would be hypersensitive to the accuracy of the computation. With A being free to 



choose, the procedure is obviously not unique. As we have remarked earlier, no unique 
continuation should be expected from a finite set of numbers, P n . However, there are 
guidelines for an optimal choice of A. 

We have mentioned in connection with (|19"D that the domain of continuation of 
q is limited by the smallest value of kj, which is k(l — A/2). Thus to increase that 
domain we would want to have a small value for A. However, with a small range 
of Xj and kj it will be necessary to have highly accurate expansion coefficients a,j in 
( |TT|) to well represent P n . So a larger value of A is preferred. This point needs to be 
demonstrated quantitatively. A way to do this is to calculate /(l) and examine its 
dependence on A. To that end we consider specific examples in the following. 

Consider a sample P n given by 

pO) = ( n + 0.3) 3 e-°- 5n /Z, n=Q,---,N (28) 

where N = 30 and Z is the normalization factor so that En=o^n = 1- For this 
V^ we have x = (n) = 7.6966 and k = 7.2193. For A chosen to be in the range 
0.1 < A < 0.6, we follow the procedure described in Sec. |3| and calculate /(l). The 
result is shown in Fig. 3(a) as a function of A. Evidently, there are fluctuations at 
small values of A, but quite stable for A > 0.25. However, when the vertical scale 
is expanded as shown in Fig. 3(b), we see that small fluctuations at a level of < 
0.3% are still present until A > 0.4, where /(l) = 7.7016. A discrepancy between 
/(l) and x is anticipated because of the finiteness of N but at 0.06% level it is 
unimportant. What is important is that A should not be too small. This example 
demonstrates the limitation of the method due to (|T9|) , if one attempts to continue 
q to large negative values. From Fig. 3 an appropriate value for A can be set at 
0.5, for which mm{kj} = 0.75 k = 5.4145. Thus F(q) can be continued to q ~ —5.4 
before encountering divergence. In practice this range of negative q is quite enough 
to exhibit the low n behavior of P n . Increasing A would improve the stability of the 
solution, but decrease the range of continued q. The choice of A = 0.5 seems like a 
good compromise between the two opposite preferences. 

Consider next another example where (n) is much smaller, corresponding to the 
situation where the phase-space cell size 5 is small. Assume 

V^ = (n + I) ' 5 e~ n /Z, N = 20 , (29) 

for which x = 0.8352 and k = 1.395. Fig. 4 shows the result of /(l) vs A, which is 
rather free of fluctuation for all A except near 0.1. The value of /(l) = 0.8352 is equal 
to x to 5 significant figures. Choosing A = 0.5 gives mm{kj} = 1.046 which does not 
allow q to go much beyond — 1. Decreasing A would not improve the situation due 
to the limitation of small k for V^ . Since A should not be changed in the analysis 
of data with varying 5, we suggest that A be fixed at 0.5. 



5 Continued Factorial Moments 

The continuation procedure having been completely specified in Sees. [| and f|, we 
can now proceed to the study of the normalized factorial moments F(q). We continue 
to use the three distributions VJfi, i = 0,1,2, as our sample inputs for P n . Unless 
otherwise stated, A = 0.5 is used. 

For the Poisson distribution V^\ let us consider the same two cases: (a) (n) = 6 
and (b) (n) = 1, already examined in Sec. ^ where the simple continuation scheme 
n\ — > T(n + 1) was used . Now, we assume N = 30 for (a) and N = 10 for (b) as 
the upper limits of n in the experimental P n . UN — oo, V^ would give F 2 = 1, 
so according to ( pip k would be infinite. For the finite N chosen for the two cases, 
F 2 are still very close to 1, so k would be extremely large. For our calculation it is 
sufficient to set k = 10 4 . Using our procedure of continuation the results are shown in 
Fig. 5(a) and (b), for the two cases (a) and (b), respectively. Note the high resolution 
of the vertical scale. For (n) = 6 in case (a) F(q) is essentially 1 for —10 < q < 30. 
In case (b) where (n) = 1, F(q) is almost 1 for —10 < q < 20 except near the edges 
of that range. This case is more difficult to continue accurately because there are 
fewer values of nonvanishing V^ . Theoretically, V^ is nonzero for any finite n, but 
we cut off at iV = 10 to simulate a realistic situation where (n) is only 1. With only 
11 values of a,j the continuation to large \q\ cannot be expected to have extremely 
high accuracy. That is why F(q) deviates from 1 near the two ends in Fig. 5(b). 
Nevertheless, the deviation is only of order 0.1%. Upon comparing Fig. 5 to Figs. 1 
and 2, the advantage of this method over the simple scheme of Sec. [| is self-evident. 

For V^ given in (|28|) with N = 30 the result of our calculation for F(q) is shown 
in Fig. 6 (a) and (b). The dependence on q is evidently very smooth. It grows rapidly 
at negative q, even though the first singularity is located at q < —5. We have also 
calculated F(q) for A = 0.3, the result of which is plotted in dashed lines, lying very 
close to the solid lines for A = 0.5. For q > the two cases are indistinguishable. For 
q < — 2 the difference is actually not negligible in absolute value but because of the 
rapid rise of F(q) it is not significant in terms of percentage discrepancy. In Fig. 6(b) 
we see that for the range of q shown the difference between the two cases is totally 
insignificant. This is the most important range for the continuous q problem, and we 
have found a reliable continuation of F(q). 

It should be pointed out that F(0) = 1 is not accidental. From (21) and (22) we 
see that F(0) = ^jdj, which is 1 by virtue of the normalization of P n and P^ B {j) in 
(11). Of course, /o is also 1 if the continuation scheme of (8) is followed. 

Finally, let us come to the third example where V^ is as given in (^9|) with 
N = 20. Now, the values of x and k are small. The calculated result for F(q) is 
shown in Fig. 7 (a) and (b). There is a fast rise at large q because of the smaller x 
(compared to the case above). The continuation to negative q encounters irregularity 
due to the small value of k. For q > —0.5, there is essentially no difference between 
the use of A = 0.5 (solid line) and A = 0.3 (dashed line). Only the solid line is plotted 



in Fig. 7(a); both are plotted in Fig. 7(b). The difference between the two A cases 
becomes noticable and quantitatively significant only for q < —0.5, a region very close 
to the singularities. The reliability of our continuation should therefore be restricted 
to the domain to the right of the sudden downturn of F(q) around q = —0.5. In that 
domain our result is smooth and insensitive to A. 



6 Multifractal Analysis 

With F(q) continuable to noninteger q, it is now possible to consider multifractal 
analysis, assuming that F(q) has a power- law dependence on the resolution scale S 
for a range of q covering both positive and negative values. Such an analysis was 
suggested previously using G moments, which are defined for all q by 

<?(?)= E(£Y (30) 

where n« is the multiplicity in bin i, nt = J2i n i, an d the sum in % is over all nonempty 
bins. G(q) shows scaling behavior 

G{q) ~ S T{q) (31) 

for q > 1 in both experimental data and model simulation, when 6(rii — q) is included 
in the summand in ( |50| ) M. However, for q < 1 (|3T| ) is not valid for any extended 
range of S, so multifractal analysis cannot be made for that range of q. The problem is 
rooted in the empty-bin effect and the fact that G(q) contains statistical fluctuations. 
Since F(q) is now defined for noninteger q, its scaling behavior 

F(q) oc 5-^ q) , (32) 



especially for negative q, should be checked for a variety of existing data. If ( |3"2"D is 
valid for a range of q around 1, multifractal analysis can then proceed without the 
necessity of subtracting out the statistical component, as was done for G(q). 

We can relate G(q) and F(q), if we assume that ( j30|) is defined for the dynamical 
distribution only without the statistical fluctuation. The sum over all bins in fl3"UD 
can then be related to averaging over the dynamical distribution in (|5]). Using (^|), 
(§JJ) and d32]), we get 



r(q) =q-l- V {q) , (33) 

where the —1 comes from the fact that YU by itself gives the total number of bins, 
which varies as £ _1 . Multifractal spectrum is then obtained by the Legendre transform 

f(a) =qa- r{q) (34) 

10 



with 

a = dr(q)/dq . (35) 

Thus the verification of the power law (^) and the capability of calculating a by 
differentiation with respect to q, which we now have, make possible the presentation of 
the scaling properties of dynamical fluctuations in terms of the multifractal spectrum 

/(«)■ 

It is not our purpose here to examine specific dynamical models or experimental 
data and to extract their multiplicity fluctuation behaviors. However, we can demon- 
strate the nature of f(a) if we assume that the scaling behavior (^) is true for a set 
of P n (S) for a range of 5. Let us further assume that among those P n (S), a specific 
one at some S is exactly V£' given in (p8|) . Then we have 

<p(q)=c\RF(q) + c 1 (q) , (36) 

where c = (— ln<5o) -1 and C\(q) is some function of q independent of 5q arising from 
the proportionality factor in (|3"2"|). Scaling behavior means that <p(q) is unchanged, 
as 5 is varied from do- The major part of the q dependence of (fi(q) derives from 
that of F(q), which we know from Fig. 6. Apart from the unknown constant c and 
the unknown function c±(q) in this example, we can determine f{o) from F(q) by 
varying q parametrically For illustrative purpose, let us assume that c\{q) = so 
that using fl33"D-(|36]) and the result of our calculated F(q) pertaining to V^\ we can 
determine f(a). In Fig. 8 we show the result for four possible values of c. We stress 
that in a model or data analysis c is not a variable; (p(q) is determined from the 
log- log plots when there is scaling. Fig. 8 merely illustrates the possible form of f(a), 
if if(q) happens to coincide with the result of analyzing a particular P n = V£' with 
a specific c in (|36"D and with Ci(g) = 0. The dashed line indicates where f(a) = a, 
and is tangent to each of the f(a) curves at q — 1. The range of q covered by /(a) in 
Fig. 8 is, depending on c, roughly — 2 < q < 4, with the q = point always occurring 
at the peak of the f(a) curve. The multifractal dimension D q is || |16| 



D q = r{q)/{q-l) , (37) 

which is related to a by 

D = f(a ) , D 1 = a 1 , (38) 

where a q is the value of a at q. Thus «o is where f(a) is maximum, and a\ is where 
j{pt\) = oti. The multifractal spectrum f(a) is the most elegant way of displaying 
the scaling properties of dynamical fluctuation. Reduced to the bare minimum, two 
parameters can be used to characterize f(a), viz. «o and a\, the location of the peak 
and a measure of the width, respectively. 



In many multiparticle production processes the scaling behavior fl32|) is not valid 



over an extended range of 5. The multifractal analysis described above cannot be 
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applied then. However, it has been found phenomenologically [ |IT| , IS, [19| as well as 
theoretically [13|, 20], not only in hadronic and nuclear collisions, but also in quantum 



optics |1J|, that F q satisfies a different scaling law 

F q <xF£" . (39) 

Let us assume that this behavior can be established for continuous q so that the 
function f3(q) can be determined for a range of q values both positive and negative. 
Then it is possible to define formally another spectrum, call it g(oe), in exact analogy 
to (|3~3"1)-([35|), but without the geometrical implication of multifractality Thus we 
define 

a(q) = q-l-(3(q) (40) 

g (a) = qa- a(q) (41) 

a = da(q)/dq . (42) 

The only input into this scheme of description is f3(q), which is 

(3(q) =lnF(g)/lnF(2) + b(q) , (43) 

where b(q) is the log of the proportionality factor in (|39|). Mathematically, a(q) and 
g(a) correspond to r(q) and f(a) if we set c = l/hiF(2), but physically (|32|) need 
not be true, rendering ip(q) meaningless, while (|39|) can well be true (no exception 
having been found so far). Ref. |13|] gives an explicit example of (|39|) being valid for 
a problem that does not have ([32]) . 

In the example where V£' is considered, let us assume that it belongs to the 
type of physical problems for which (|39"D is valid. Then the function F(q) obtained 
is sufficient to determine g(a), assuming b(q) = 0. The result is shown in Fig. 9(a). 
Since F(2) = 1.14, g(a) corresponds to f(a) in Fig. 8 with c = 7.64. In Fig. 9(a) the 
peak occurs at ao = 1.45 while the tangent point is at a,\ = 0.52. If, on the other 
hand, Vffl is used as an example that has a scaling behavior (|39"D, totally unrelated 
to V^\ then the corresponding spectrum g{a) is as shown in Fig. 9(b). Note that in 
this case the maximum a is 2.15 corresponding to q = —0.24. Our continuation to 
lower value of q has led to a sudden downturn of F(q) around q = —0.4, exhibited in 
Fig. 7(b). That causes a drastic change in the derivative in ([42]) at around q = —0.24 
which in turn gives rise to an irregular behavior in g(a) at a = 2.15. Thus the nature 
of V^ prevents the use of F(q) for q < —0.24, thereby setting an upper bound to 
how far to the right the spectrum g(a) can be developed. In that figure a = 1.76 
and «i = 0.41, quite different from the corresponding values in Fig. 9(a). These two 
figures are sufficient to indicate that V^ and V^ cannot belong to the same class 
of F-scaling factorial moments. 

In summary, g(a) is a representation of the characteristics of the F-scaling be- 
havior, (|39"1), of F(q), and may be a generally useful description of all multiparticle 
production processes. 
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7 Conclusion 

We have presented a way to determine F(q) for continuous q such that it is 1 for all q if 
the input distribution is Poissonian, i.e. the statistical fluctuation is filtered out. The 
range of q into which F(q) can be continued depends on the nature of P n . Generally 
speaking, the low n part of P n is characterized by the negative q region of F(q). Thus 
the study of the scaling behavior of multiplicity fluctuations can now be extended to 
dips, gaps and voids. All existing data that have been put to intermittency analysis 
at positive integer q should be reanalyzed for continuous q. Similarly, such reanalysis 
should be done for all models and MC codes. Thus the confrontation between theory 
and experiment can now be extended to a significant portion of the real line of q, as 
compared to the previous situation where it has been done for only a few isolated 
points at the integer values. For comparison, the study of Bose-Einstein correlation 
is focused on only one point: q = 2. 

If the dynamics of particle production is self-similar so that F(q) exhibits a power- 
law dependence on the resolution scale 5, then the intermittency index <p(q) can be 
determined as a continuous function of q. It then follows that the multifractal spec- 
trum f(a) can be derived without any ambiguity or need for correction to eliminate 
statistical contamination. On the other hand, if there is no power-law dependence on 
5, but there is F-scaling, i.e., F g ex F^ , which is more commonly observed, then the 
knowledge of f3(q) is sufficient to determine the spectrum g(a) that gives an excellent 
description of the self-similar behavior of the dynamics of fluctuations. 

For the purpose of providing a convenient outline of the procedure to determine 
F(q), f(oc) and g(a), we summarize here the steps needed to do the analysis: 

1. Starting with the input P n , n — 0, • • •, N, determine x and k from ([23|) and 
d|), and then Xj and kj from (|5|)-(|27|) with A = 0.5. 



2. Using (0), set up iV + 1 linear algebraic equations, flTT|), with P„ (j) as the 
matrix and P n as the input vector. Solve for a,j, j = 0, • • •, N. [The use of 
MATHEMATICA has been found to be convenient.] 

3. Use (^) to calculate F(q) for a range of real q. 

4. If the input P n (5) is known for a range of S, determine F(q) as a function of 5 
and examine the validities of (|32|) and (|39"D over that range of S. 



5. If (p2\) is valid for a subrange of 5, then from ip(q) and (|33| ) - (|35|) determine f(a) 

6. If ([52|) is invalid but (|59| ) is valid, then from j3(q) and (|4TiD-([42j) determine g(a) 

7. Compare theory and experiment at any of the three levels: F(q), /(a) or g(a) 
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Figure Captions 

Fig. 1 Normalized factorial moments F q in the simple continuation procedure using 
(|8|) with Poisson distribution (|9|) as input and with (n) = 6. (a) —10 < q < 25; 
(b) 25 < q < 30. 

Fig. 2 Same as Fig. 1 but for (n) = 1. 

Fig. 3 Dependence of factorial moment f(q — 1) on A when the input distribution is 
V^'. The resolution in (b) is higher than that in (a). 

Fig. 4 Dependence of f(q = 1) on A when the input distribution is V^\ 

Fig. 5 Normalized factorial moments F(q) with high vertical resolution for Poisson 
distribution with (a) (n) = 6, and (b) (n) = 1. 

Fig. 6 F(q) with V^ as input distribution. Solid line is for A = 0.5 dashed line for 
A = 0.3. (a) -3 < q < 10; (b) -1.5 < q < 2.5. 

Fig. 7 F(q) with V^ as input distribution. Solid line is for A = 0.5 dashed line for 
A = 0.3. (a) -1 < q < 6; (b) -1 < q < 2. 

Fig. 8 Multifractal spectra f(a) for four values of c (see text), if the input V^ belongs 
to a class of scaling distributions. Dashed line is for f(a) = a. 

Fig. 9 Spectrum function g(a) if (a) V^ and (b) V^ belong to separate classes of 
self-similar distributions satisfying F-scaling (3.9). 
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